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The many-body theory of degenerate systems is described in detail. More generally, this theory 
, applies to systems with an initial state that cannot be described by a single Slater determinant. 

The double-source (or closed-time-path) formalism of nonequilibrium quantum field theory is used 
to derive an expression for the average value of a product of interacting fields when the initial state 
is not the vacuum or a Slater determinant. Quantum group techniques are applied to derive the 
. hierarchy of unconnected Green functions and the hierarchy of connected ones. 
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j— — ]■ Degenerate systems are plenty (all systems containing an odd number of electrons, by Kramers theorem 1 ^) and 
degeneracy plays an important role in numerous interesting physical effects (e.g. magnetism or superconductivity). 
I ■ Thus, it seems relevant to develop calculation methods for degenerate systems. The density functional theory of 
£j \ degenerate systems is a subject of continued interes t 3 ! 4 ! 5 ! 6 ! 7 . Green functions can be useful because they provide 
an exact expression for the exchange and correlation potential^ and for their ability to calculate excitations, e.g. 
through the GW approximation^ or Bethe-Salpeter equatio n 11 ! 12 . Here, we study the Green functions of degenerate 
systems. 

The quantum field theory of degenerate systems has been investigated since the sixties by two methods: either 
'-^J \ by calculating the ^-matrix elements between different "in" and "out" statesi 3 * 1 ^ 5 ^*^ or by assuming that the 
interacting ground state is a pure state evolving from a non-interacting pure statei&iiiSSiSii 2 ^. However, these works 
treat the electron-electron interaction as a perturbation, and we know that this is a rather crude approximation. So 
i we need a non-perturbative approach to the Green functions of degenerate systems. The simplest and most common 

non-perturbative method is self-consistency. Therefore, we shall develop in this paper a self-consistent calculation the 
OnJ , Green functions of a degenerate system. 

The first problem that we meet with such a program is the fact that we cannot describe the system with a 
wavefunction. For example, the electronic configuration of a boron atom in the ground state is ls 2 2s 2 2p 1 . The 
Hamiltonian is invariant by rotation and, if the spin-orbit coupling is neglected, the six pure states \pi,s) (with 
i = x, y, z and s = ±1/2) are degenerate and are eigenstates of L 2 and s 2 . However, none of these pure states gives a 
spherically symmetric electron density. More generally, if the ground state of a quantum system is a pure state with 
£f) ' angular momentum L > 1, the charge density derived from this state is not spherically symmetric^ 3 .. Therefore, the 
self-consistent pure state (with L > 1) of a spherical Hamiltonian breaks the spherical symmetry of the problem. A 
related results was proved by Bach et ali^: the solution of the unrestricted Hartree-Fock equations does not contain 
unfilled shells. Therefore, for the boron atom, the 2p shell is deformed to lift its degeneracy. 

To cure this defect, we have to assume that the boron atom is in the mixed state described by the density matrix 
p = (1/6) \pi, s)(pi,s\, which preserves the rotational symmetry of the system. 

Therefore, for degenerate systems, we need to calculate an evolution starting not from a single pure state (usually 
called the vacuum |0)), but from a density matrix p. The best tool to do so is nonequilibrium quantum field theory, 
O ■ as created by Schwinger, Kadanoff, Baym and Keldysk 25 - 26 - 27 . In particular, the closed time-path method will enable 
us to express the various Green functions we need as functional derivatives with respect to external sources. 
• <-h , We describe now the main result of the paper. The differential form of the Green function hierarchyS 6 ^ is 

X : 

(^ + ^)G(l,l') = 5(1-1')-* /«(r 1 -r 2 )G 2 (l,2;l',2+)dr 2 . 
For non-degenerate systems that can be represented by a Slater determinant, the integral form of this equation is 

G(l,l') = GoCl.l')-* / G (l,3Mr3-r 2 )G 2 (3,2;l',2+)dr 2 dr 3 , 
where Go(l, 1') is the Green function for the free Schrodinger equation 

(^ + ^o(MO = ^(i-iO. 
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The unperturbed Green function G°(x,y) is given by the following expression (see2S p. 124) 

G (x,x') = -i6{t-t') e- ,f - (H '»%(r)%(r') 

e„>ep 

+i9(t'-t) ]T e-^-^Uni^Un^'), (1) 

with x = (t, r), x' — (t' , r') and the orbital u n (r) is the solution of the unperturbed Schrodinger equation for energy e n . 
The Fermi energy tp is chosen so that the total charge — i J tr(Go(t, r, t, r))dr is equal to the number N of electrons 
in the system. In this independent-particle picture, the ground state is degenerate if the Fermi level is degenerate and 
not completely filled. Therefore, the definition Q must be modified because it assumes that the Fermi level is full. 
We shall sec that the definition of Go for a degenerate system is non trivial. 

Moreover, for degenerate systems, the relation between the differential and the integral Kadanoff-Baym equations 
is modified, because the solutions of the free Schrodinger equation intervene. In fact, the full hierarchy of Green 
functions is changed. 

The correct hierarchy of Green function is important because it is the basis of the GW approximation 1 ^. Thus, the 
GW approximation must be adapted to degenerate systems. A similar modification is required for the Bethe-Salpeter 
equation. 

In this paper, we give the proper expression for Go and the integral form of the Kadanoff-Baym equation for a 
general density matrix. Compared to previous results, the present equations have two advantages: they are adapted 
to a self-consistent treatment and they do not break the symmetry of the problem. 

Although the quantum field theory of degenerate systems seems to be a rather natural problem, it was not solved 
before because it poses technical difficulties that can hardly be overcome with the standard many-body techniques. 
Our main tool here will be the quantum group (or Hopf algebra) approach to quantum field theory, developed in 30 31 . 

In this paper, we give a self-contained presentation of the calculation of the expectation values of products of quan- 
tum field operators in the interaction representation. Then we compute interacting Green functions using functional 
derivatives with respect to external sources. The Hopf algebra of derivations is then introduced and used to derive 
the hierarchy of Green functions for systems with degenerate initial states, or more generally for systems with initial 
correlation. Explicit hierarchies are obtained for unconnected and connected Green functions. In a forthcoming pub- 
lication, the special case of a single electron in a system with closed shells and a two-fold degenerate orbital will be 
calculated in detail. 



II. EVOLUTION OF EXPECTATIONS VALUES 

We saw in the introduction that a self-consistent calculation of degenerate systems requires the use of density 
matrices. In this section we investigate how the unperturbed density matrix evolves with time under perturbation. 
As a first step, we calculate the evolution of an unperturbed wavefunction, then we extend this to the evolution 
of a density matrix, and we use this result to calculate the evolution of an expectation value. The calculation of 
transition amplitudes in quantum field theory is not completely standard, so we give here a detailed derivation. As 
an application, we obtain a formula for the Green function of a degenerate system. 

A. Evolution of wavefunctions 

We start from a time-independent free Hamiltonian Hq — J ?/'5(r)/io(r)V's(r)dr, where ho(r) is a one-particle Hamil- 
tonian and ips( r ) is the field operator in the Schrodinger picture. A convient form of ips( r ) is , 0s( r ) = ^2 n u n( r )b n , 
where u n (r) is an eigenstate of Hq: hou n = e n u n , and b n is the annihilation operator for the one-particle state u n (r). 
We first look for the solutions of the Schrodinger equation 

= ffo|#°(t)>. 

As usually, we isolate the time dependence by putting = e~ l£ °*|$°) so that if |$°) = ^"l^n)- We assume 

that the |$° } provide a complete set of states. The matrix elements of the operator As(t) in the Schrodinger picture 
is ($^(i)|As(i)|$° (£)). We go now to the Heisenberg picture by 

= <eie^s(t)e-^|$°> = (<S>° m \A(t)\*° n ), 
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where A(t) = e lHot As(t)e lHot is the operator As(t) in the Heisenberg picture. In particular e lH ° t Hoe lHot — 
H , so that H is the same in both picture. The field operator in the Heisenberg picture is (see^ p. 146) ip(x) — 
e iHot ip s (r)e- iHot = E„ M r ) e ~ ie "* & n> with x = (t,r). 

We are interested in the interacting theory, so we add a possibly time-dependent interaction term to the Hamiltonian. 
This gives us Hs(t) = Ho + Hcj^ty) in the Schrodinger picture and H(t) = Ho + H lnt (t) in the Heisenberg picture, 
with H iDt (t) = e iffo 'i3^ nt (t)e" lffot . In practice, H^Jt) is a polynomial in ip s (r) and V>s(r), and H int (t) is the same 
polynomial where V's( r ) is replaced by tp(t,r) and V>s( r ) is replaced by ip(t,r). 

We look for solutions of the Schrodinger equation 

ij^W) = H s (t)\<t> s n (t)). 

We go to the Heisenberg reprensentation with respect to Ho (which is called the interaction picture) and we define 
\<f> n (t)) =e iHat \$%(t)). Therefore, 

i^|*»(*)> = e iH ° t (-Ho + H s (tm S n (t))=H i » t (t)\<t> n (t)). 

To solve this problem, we look for an operator V(t) such that | $„(£)) = V(i)|$°). The Schrodinger equation becomes 

= H^(t)V(t)\<S>° n ). 

This must be true for the complete set of thus 

i^V{t) = H^(t)V(t). (2) 

B. Calculation of V(t) 

To solve equation 0, we put U{t,t') = V{t)V~ l {t'). Therefore, U{t,t) = 1 and 

i§-U(t,t>) = H^(t)U(t,t>). 

We are going to prove some properties of U(t, t'). We first prove the group property U (t, t')U(t' , t") — U (t, t"). From 
the fact that V(t)V -x (t) = 1 we deduce 

8 



and 



Thus 



l di V {t) = ~ V ~ (t)H W ' 



i^U{t,t') = -U{t,t')H™\t% 



i-^U(t,t')U(t',t) = U(t,t')(-H i « t (t') + H™ t (t'))U(t>,t) = 0. 

Hence, the product U(t,t')U(t' ,t) is independent of t'. To find its value, we put t' — t, so that U(t,t')U(t' ,t") 
U(t,t)U(t,t") = U(t,t"). Then we show that U(t,t') is unitary. We take the adjoint of equation 0: 

-*J^t(t) = W^ir^Ct), 

because H lnt (t) is Hermitian. This implies 

i^ f (*»0 = K- lt (* , )*^ t (t) = -^t(t,f)f^( t ). 
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Therefore, W(t,t') = U(t',t) because both operators satisfy the same equation and the same boundary condition 
U(t,t) = W(t,t) = 1. But the group property leads to U(t,t')U(t' ,t) = U(t,t) = 1, so that W{t,t') = U{t',t) = 
[/^(MOi and U(t,1?) is unitary. 

The construction of U(t,t') is standard (see, e.gi 32 i 33 ) and yields 

U(t,t') = Texp(-i^" H int (T)dr). (3) 

Here, T is the time-ordering operator that orders its arguments by decreasing time from left to right. For example 
T(A(t)B(t')) is A(t)B(t') if t > t' and is B(t')A(t) if t' > t. An important property of the time-ordering operator is 
that its arguments commute. For instance, it can be checked from the definition that T(A(t)B(t')) = T(B(t')A(t)). 
To complete the picture, we use the adiabatic hypothesis which states that 

lim |$„(t)> = |$ ), 

t — > — oo 

so that 

lim V(t) = 1, 

t — > — oo 

and V(t) — U(t, — oo). Thus, V(t) is unitary. This has two important consequences: (i) the states | $„(£)) are complete 
at all times: 

£l<M<)K*«(t)l = v{t){Y,\*° n )K\)v\t) = v{t)v\t) = i, 

n n 

and (ii) the scalar products are conserved: ($ m (i)|$ n (i)) = ($^|$°}. 

To complete this section, we define the notion of anti-time-ordering operator. For any X which can be written as 

a product of field operators, the anti-time-ordering of X is defined as T*(X) — (T(X^)y (see 3 ^). Notice that T* is 
linear and its arguments commute. To understand the physical meaning of T* , we take an example. If t > t' we have 
T*(A(t)B(t')) = (T(Bt(f)At(t))) t = (A^(t)B^(t')) f = B(t')A(t). Analogously T*(A(t)B(t')) = A(t)B{t') if t < t' . 
In other words, T* orders its arguments so that the operators are on the right when they occur later. This is true 
for any number of arguments, and T* orders its arguments in the reverse order with respect to T. This is why T* is 
called the anti-time-ordering operator. The main example is 

T*(V>(a#0/)) = %° - x )i>(x)i>(y) - 8(x° - y°)j(y)^(x). 

The most important application of the anti-time-ordering operator is the calculation of U^(t,t'). 

U\t,t') = (Texp(-^ i i/ int (r)dr)) 

= T*exp(iJ iJ int (r)dr), (4) 

because 7J mt (r) is Hermitian. 



C. Evolution of density matrices 

If !$"(£)) are solutions of the Schrodinger equation for Ho, a density matrix p%(t) in the Schrodinger picture has 
the following general form 

ran 

where p nm is a Hermitian matrix with non-negative eigenvalues such that p nn = 1. For later convenience, we do 
not require p nm to be a diagonal matrix. From the Schrodinger equation, we see that the density matrix satisfies the 
equation 



— tt— - -t[H ,p s (t)\. 
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As for the wavef unctions, we define the density matrix in the Heisenberg representation p by p = e %Hot p%{t)e~ ot = 
Ylmn Pnrn\<&n){® m \- Notice that p does not depend on time. 

In the interacting case, we look for a density matrix ps(t) in the Schrodinger picture that we write 

P S (t) = £p™l*n(*)><*m(*)l- 
mn 

It satisfies the equation 

dps(t) 
dt = -i[H,Ps(t)\- 

We go to the interaction picture by defining pi(t) — e lHat ps(t)e~ lHot , which satisfies the equation 
Now it is easy to see that the density matrix 

Pl(t) = Y,P^V(ml)(K l \VHt)=Y,Pnrn\^n(t))(^ m (t)\, 

mn mn 

satisfies the above equation. In other words, the density matrix pi{t) describes the interacting system and it can be 
considered as the interacting density matrix that evolved from the non-interacting density matrix p because of the 
interactions. 

D. Evolution of expectation values 

The value of the observable A(t) (in the interaction picture) for a system in a mixed state described by the density 
matrix pi is (see^ p. 314) 

(A(t)) = tx(piA(t)) =J2Prnn(®m(t)\A(t)\$ n (t)) 

III 1 1 

= ^p^($°JV*(t)A(t)V(t)\*° n ) = tr{pV\t)A(t)V(t)) 

in ii 

= tv(pU(-oo,t)A(t)U(t,-oo)). (5) 

The group property of U{t, t') enables us to derive 

(A(t)) = tr(pU(-oo,t)U(t,+(x)U(+(x>,t)A(t)U(t,-oo)) 
= tx(pU(-oo,+oo)U(+oo,t)A(t)U(t, -oo)), 

= tr(pS- 1 T(A(t)e- iA! * t )), (6) 

where the interacting action is (up to a sign) _4 lnt = i/ lnt (r)dr and where the S-matrix is defined by S = 
U(+oo, — oo) = T(c~ lA ' nt ). The last line of © was derived as follows. By equation (J3J) 



U(+oo,t)A(t)U(t,-oa) = T(exp(-i J H mt (t)<1t)) A(t) 

T(exp(-i I ff int (r)dr)). 

J — oo 

In that expression, the operators are on the left when their time arguments are larger. Thus, they are time ordered 
and we can rewrite this 



U(+oo,t)A(t)U(t,-oo) = T^expH H iut (T)dT)A(t) 

expH / H int (T)dT)). 
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The arguments of the time-ordering operator commute, thus 

U(+oo,t)A(t)U(t,-oo) = T^A{t)exp{-i J # int (r)dT) 

exp(-i / H mt (T)dT)) 

J — OO 

= r(A(t)exp(-< f H hlt (T)dT) 
= T(A(t)c- lAint ). 

To obtain equation ©, we inserted 1 = U(t,+oo)U(+oo,t) before A(t) in equation (J5J. Of course, we can also 
insert 1 = U(t, +oo)U(+oo,t) after A(t) in equation J5j. This gives us the alternative formula 

(A(t)) = tr(pT*(A(ty Aiat )S). (7) 
E. Correlation functions 

Finally, we shall have to determine the correlation function between an observable A(t) at time t and an observable 
B(t') at time t' . To do this, we must determine which picture must be used to describe the observables at two different 
times. It turns out that the Heisenberg picture does the job. There are three reasons for this: (i) the equation for the 
observables in the Heisenberg picture are similar to the equations for the corresponding classical observables (see^ 
p. 316), (ii) the correlation functions of observables calculated in the Heisenberg picture agree with the experimental 
measurement of these observables (seeS, p. 655), (iii) the quantum description of photodctcctors shows that they 
measure the correlation functions of the photon field in the Heisenberg picture (see2£, chapter 14). 

The relation between the observables in the Schrodinger and Heisenberg pictures is given by Ajj(t) — 
Vg(t)As(t)Vs(t) (see^, p. 143), where Vs satisfies the Schrodinger equation for the full Hamiltonian Hs(t): 

^ = -iHsWs®. 

The standard boundary condition is Vs(0) — 1 and the solution of this equation is Vs(t) = c~ lHot U(t, 0). The boundary 
condition means that the Heisenberg and Schrodinger pictures coincide at t = 0. Therefore, the time-independent 
density matrix of the Heisenberg picture is equal to the Schrodinger density matrix at t = 0, i.e. 



ph = Mo) = £aw|*2(o))«(o)|. 

ran 

The correlation function for the two variables A(t) and B{t f ) is now 
(A(t)B(t')) = tr(p H A„(t)B H (t')) 



mi i 



= Y / Pnm« l \U(-^,0)A H (t)B H (t , )U(0,-oo)\'i> n ) 

mn 

= J2 Pn m ($° m \U(-°o, 0)17(0, t)e iHot A s (t)e- iHot U(t, 0) 

mn 

U(0,t')e lHot 'Bs{t')e~ mot 'u(t',0)U(0,-oo)\^° n ) 
= tr(p[/(-oo, t)A(t)U(t, t')B{t')U{t', -co)) . 

As in the previous subsection, the group property of the evolution operators U(t, t') enables us to rewrite three kinds 
of correlation functions, for the operator product of fields, the time-ordered product of fields and the anti-time-ordered 
product of fields. 

(A(t)B(t')) = tr(pT*(A(t)e^' nt )T(i?(t')e- l ' 4int )), 
(T(A(t)B(t>))) = tr(p5- 1 T(^(t)fl(t')e- wtot )) ) 
(T*(A(t)B(t'))) = tr(pT*(A(t)B(t')e lAint )S). 
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III. FUNCTIONAL DERIVATIVE APPROACH 
A. Functional derivatives of the S-matrix 

The use of functional derivatives in quantum field theory was advocated by Schwinger 36 . The S-matrix for a 
nonrelativistic systems of electrons with Coulomb interaction is given by 

S = [/(+oo,-oo) =T(e- lAiat ). 
In solid-state physics, we usually consider the free and interaction Hamiltonians ( 37 p. 44) 

2 

H ° = E / $s(t,r)(-£- + U N (r))Mt,r)dr, 

S—l 

# int W = \Y, /^(^r)^(t,r')V;(r-r')^(t,r')^(t,r)drdr', 

s,s' 

where Un(y) describes the interaction with the nuclei and V e (r) = e 2 /(47reo|r|) the electron-electron interaction. We 
define now an S-matrix which depends on two external fcrmion sources rj(x) and fj(x) as 



5(77,77) = T exp ( — iA mt + i J r)(x)ip(x)dx + i J ip(x)r)(x)dx) . 

For a nonrelativistic fermion, ip(x) and ip(x) are two-component vectors. Thus, the sources are also two-component 
vectors and 

2 2 
fj(x)ip(x) = ^2fj s (x)^p s (x), $(x)i](x) =^2?/j s (x)r ls (x). 

s=l s=l 

The functional derivative with respect to the fermion source 77(2) satisfies 

6 S 

My) = 8(x-y), ——f](y) = 0, 



8rj{x) ' 5rj{x) 



5r)(x) 5r)(x) 5r)(x) 

In this equation, we assumed that u is the product of a certain number of fermion fields or sources, and this number 
is denoted by \u\. Similar relations are satisfied by the functional derivative with respect to fj(x). Equation JSJ is 
known as Leibniz' rule. 

The sources 77 and 77 anticommute, so the functional derivatives anticommutc: 



5r]{x)8r](y) 8r)(y)5r)(x) 

To see how functional derivatives act with respect to the time-ordering operator, we first notice that the sources 
can be taken out of the time-ordering operator. For example, if x° > y° 

T(fj(x)i/j(x)4>(y)ri(y)) = fj(x)ip(x)ij(y)r](y) = f](x)f](y)ip(x)ip(y) 
= f){x)r)(y)T(Tl)(x)$(y))i 

if x° < y° 

T(fj(x)rp(x)ip(y)r)(y)) = tp(y)ri(y)fj(x)ip(x) = -f){x) i q{y)i}{y)tp{x) 
= v(x)v(y)T(^(x)i[>(y)). 

Thus, the functional derivative with respect to 77(2:) or fj(x) commutes with the time-ordering operator. In particular 

^'^-l^o = iT(^(x)e- lAiat ) =iU(+oa,t)il>(x)U(t,-oo), 



Srj(x) 

''^'fi=,=o = -iT(i>(x)e- iA ™) = -iU(+oo,t)$(x)U(t,-oo), 



5r)(x) 
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where x = (t, r)2£ and the minus sign in the last equation comes from the fact that the functional derivative must 
jump over xj}(x) to reach f}{x) in the definition of S(fj,rj). With this definition, we can write 

(Mx)) = i^rE^-WI 5 ^^)" 1 ^^)!^)!^^. (9) 

' ^ ' ran 

In the vacuum, the density matrix is |0)(0| and 

& H (x)) = i__(0|£(0,0)- 1 £fo rj) |0)| fl=fl=0 . 



One then invokes the "stability of the vacuum"- to derive 



6 (0\S(fj,v)\0) , (10) 

which is the Gell-Mann and Low formula 39 . The denominator is a pure phase, thus the main problem is to calculate the 
numerator of equation (|10|l . A standard result of the functional derivative approac h 3 ?i 4Q is that the interacting S-matrix 
S(fj, rj) can be obtained from the non-interacting S-matrix S°(fj, rj) with S°(fj, rj) = Texp (i J fj(x)ip(x) + *ij)(x)r)(x)dx) 
by the equation 



/OO ■ C ■ C 



where x = (t, r). For a state described by a density matrix p — p„„ 1 |5 , q)($™|, the Gell-Mann and Low formula does 
not hold and we must deal with the term 5(77, ?7) — 1 in equation This is done by doubling the sources. 



B. Source doubling 

The idea of doubling the sources was proposed independently by Schwinger— and Symanzik^L^. It is a basic 
technique of nonequilibrium quantum field theory^ 3 ^^^^^^ 9 " where it is also known as the closed time-path 
Green function formalism. For equilibrium quantum field theory, Wagner showed that it can be useful to triple the 
sources^,. In equation JUJ, we have the operator product of S(fj,7j)~ 1 and S(fj,rf). We cannot obtain an operator 
product by functional derivatives, because they generate time-ordered products of operators. Therefore, we shall use 
sources to calculate S(fj,r])~ 1 and sources to calculate S(fj,rj): we define 

Z P = Y,P^o\S(fj-,rj-r 1 S(fj + ,r 1+ )\^). (11) 

mn 

Here Z p is a function of the sources fj- 1 rj-,fj + .rj + . Notice that Z p = 1 when fj- — fj + and r/_ = rj + , because 
S{fj 1 rj)~ l S{fj,rj) = 1 and trp = 1. To calculate S(fj,r])~ 1 , we recall that S is unitary, so that 

SiV-,!-)^ = S(fj-,r]-)^ 

= (Texp ( - iA int + i J fj-(x)i>(x) +$(x)r)-(x)dx)y 

= T*exp (iA mt - i I fj-(x)ip(x)+ip(x)r)-{x)dx), 



where T* is the anti-time-ordering operator first considered by DysonSLS (see alsoS p. 94), which orders operators 
according to decreasing times. For example, 

T*(i>(x)jj(y)) = %° - x a )il>{x)${y) - 6(x° - y a )${y)i>{x). 

As for S(fj+ : 77_|_), we can write 

/OO > £ • C 

-00 5r)-{x) 6r)-{x) 
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and 



where x — (t, r) and 

= r e cp(-i/,-_ WW d,-i/^_(^). 

If we put all this together, we obtain 

Z p = e- iD Z° p , (12) 

where 

i-oo <"7+(z) Sr) + {x) 6r/-{x) 6r)-(x) 

Z° p = Y,Pn m {< l \S°(^,^r 1 S°(v + ,r 1+ )\^ l ). (14) 

run 

Notice that the functional derivatives with respect to r)- (x) and fj^ (x) correspond to anti-time-ordering. 
These are the basic equations for the calculation of Z p . The next step is now the evaluation of Z®. 

IV. CALCULATION OF Z° 

In the calculation of Z p , we first write S (fj_,r]_)~ 1 S (fj + ,r) + ) in terms of normally ordered operators, then we 
calculate the trace of the normal ordered term. The use of normal order is very convenient to calculate matrix 
elements. 

A. Normal ordering 

If we call A = —if fj^(x)i/j(x) + ip(x)r]_(x)dx and B = if fj + (x)il>(x) + ip(x)r] + (x)dx, we have 
S°{fj-,V-y 1 S°(f] + ,7 1+ ) = T*(e A )T(e B ). We want to write T*(e A )T(e B ) as the product of scalar terms with the 
normally ordered exponential :e A+B :. To achieve this, we use the identity giving the time-ordered exponential in 
terms of the normally-ordered exponential: T(e B ) = e^:e B : (see eq.(4-73) p. 183 of refi^S), where 

(3 = ~ j ' fj + (x){Q\Tip(x)$(y)\0)r] + (y)dxdy. 

This identity is a generating function for Wick's theorem. The same proof leads to the corresponding identity for the 
anti-time-ordered products T*(e A ) = e a :e A :, where 

a = — J fj-(x){0\T*^p(x)^p(y)\0}r]^(y)dxdy. 

Thus, T*(e A )T(e B ) = e a+ @:e A ::e B : and it remains to normally order the operator product of :e A : and :e B :. To do 
that, we write the operator exponential in terms of a normally ordered exponential e A — e a :e A : and e B = e* 9 :e B :, 
where 

a' = -\ [ fj-{x)(0Mx),^)]\0)ri-(y)6x6y, 

& = -\ J fj+{x){0\[i;(x)M]\0)Ti+(y)dxdy. 

This identity is the generating function for Wick's theorem for operator products. To obtain this result we start 
from eq.(4-72) p. 183 of ref. 38 and we use the fact that (x), ^ (+) (y)} = (0\^(y)ip(x)\Q) and {4>^(x), ipW (y)} = 

{0\ip(y)ip(x)\0) . Thus, :e A : :e B : — e~" ~@ e A e B . To transform the product e A e B , we can employ the classical expression 
e A e B = e A+B+[A,B]/2^ valid when [A, B] commutes with A and B (eq. (4-15) p. 167 of ref. 38 ). This is the case here 
because 



[A,B] = I fj-(x){^(x)^(y)}T ]+ (y) 

-V-{x){ip(x),tp(y)}fi+(y)dxdy, 
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is not an operator but a function (i.e. {ip(x),ip(y)} = (0\{ijj(x),ip(y)}\0) . Now, we transform again the exponential 
e A+B into a normally ordered exponential by e A+B = e 7 :e A+B :, where 



Vd(x)(0\[%jj(x),^(y)]\0)rid(y)dxdy, 



with rjd = rj + — i]- and T)d = t]+ — T]-- Putting all this together, we find :e A ::e B : = e^e" 4 "*" 5 :, with 7 = —a' — j3' 
[A,B]/2 + 7', so that 



7 = / V-(x){0\ip{x)iJj(y)\0)fj + (y) 

+fj- (x) (0| i>{x)i){y) |0}?7 + (y)dxdy. 
Thus, T*(e A )T(e B ) = e Q +' 3+7 :e A + iJ :. The calculation of a + (3 + 7 gives us 

T*(e A )T(e B ) = expH / v(x)G Q (x,y) V (y)dxdy}N (fj d , Vd ). 

The two-dimensional vectors 77 and 77 are 

«*>-( SrfS ) • 

the free Green function is 

r o rr s _ ^-<(0|T(^(x)^(i/))|0) -i(0|^(tf)^(a:)|0> \ , . 

ot - ' yJ ~ V i(Q\i>(x)i,(y)\0) -i(Q\T*(i;(x)i;(y))\0) J > [W) 

and the normally ordered exponential is 

N°(f}d,r) d ) = :exp [if fj d (x)'ip(x) + $(x)i] d (x)dx]:. 
Notice that the Green function is a solution of the equation 

{id x o - h Q (x))G° (x, y) = ( 6{X V) _ (5(x °_ y) ) . 
Finally, the generating function is 

Z° = expH f fj(x)G (x,y)r](y)dxdy}tr[pN (f] d ,r]d)}. 



A similar expression is given in ref4&. 

Schwinger 25 showed that this expression can be rewritten in terms of advanced and retarded Green functions, using 
the sources fj m — (fj + + fj-)/2 and r/ m = (j] + + rj-)/2. 

Z p = cx Ph J Vd(x)G r (x,y)rj m (y) - fj m (x)G a (x,y)r) d (y) 
1 - - -i0/„ .A m ./-..\j„j..u.rsRrO/ 



with 



+-^Vd{x)G° c (x,y)r]d(y)dxdy]tT[pN (%,%)], 



G°(x,y) = 6(x -y°)(Q{ip(x),Tp(y)}\0}, 
G° a (x,y) = -6(y -x )(0\Mx),i,(y)}\0), 
G° c (x,y) = (0\mx),i>(y))\0). 
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B. Calculation of tr[pN° (fj d ,Vd)} 

The calculation of tr[pN°(fjd, rjd)} is relegated to appendix because it is rather technical. We give here the results. 
The unperturbed eigenstates of H will now be called \K) and \L) instead of |$^) and |3>°). They are defined from 
the vacuum |0) by application of creation operators \K) — b\ . . . |0) and \L) — b^ N . . . b'^ |0). Here, N is the number 
of electrons and the indices ik and jk take their values in the set of indices of the M orbitals. We assume that the 
indices are ordered: i\ < • • • < ijy and ji < • • • < Jn- If we take the example of Cr 3+ , the number of d electrons 
is N — 3 and the number of d orbitals is M — 10. We assume that the orbitals are ordered in such a way that 
the M orbitals that come into play are numbered from n = 1 to n = M . We define integrals of the product of the 
wavefunctions with external sources by a n = J fjd{x)u n (x)dx and a n = J u n (x)r)d(x)dx, where u n (x) = e _le " t it„(r) 
and u n (x) = e le ™ u n (r), with x — (t, r). Recall that fjd = f}+ — fj- and rjd = T}+ — f)—- The result can now be stated 
in its simplest form as ti(pN°(fjd,r]d)) — J2kl PlrN'j {l with 

N° KL = (K\N (fjd,V d )\L) 
M d 2 

= ^(E^;)^^-^' ( 16 ) 

n— 1 

A more explicit but more cumbersome form of this result is given in the appendix. 

It is interesting to consider the particular case of a closed shell (see equation l|56ll in the appendix). This happens 
when all orbitals are occupied, i.e. N — M. Then there is only one state, p = 1 and 



tr(piV°) = + a ih a ik ). 



k=l 



V. THE HOPF ALGEBRA OF DERIVATIONS 



The term quantum group has a broad meaning^ 4 -, ranging from general Hopf algebras to g-deformed groups. In 
this section we use the more precise term of Hopf algebra. 

We introduce now the Hopf algebra of functional derivations P, which plays a vital role in this paper. In particular, 
the calculation of tr(pN°(fjd, rjd)) and the resummation leading to the hierarchy of Green functions for degenerate 
systems make essential use of the Hopf structure of V. Writing this hierarchy without Hopf-algebraic tools would be 
quite cumbersome. Since the introduction of the Hopf algebra of renormalization by Kreimer— , it has become clear 
that Hopf algebras are going to play a substantial role in quantum field theory22i2i. 

Many textbooks on Hopf algebras are now available^S but we shall use only a very limited amount of this theory. 
For the convenience of the reader, we give now a short survey of the Hopf algebra of derivations. 



A. A familiar example of coproduct 

The most unusual object of a Hopf algebra is the coproduct. To make the reader familiar with this concept, we 
present it in the case of the algebra A of differential operators with constant coefficients. We consider the coordinates 
X\, . . . ,x n of an n-dimensional space, and the differential operators P — ^2 a a a D a , where a — {a\, . . . , a n ) is a 
multi- index, a a is a complex number and D a — . . . 9™ 1 , where di denotes the partial derivative d/d Xi . It is clear 
that A is a vector space with basis D a , where a runs over all the possible multi-indices. A is also an associative 
algebra with the product induced by the product of the basis elements D a D" = D a+ " . To this algebra we add a 
unit 1 such that Dl = ID = D for any element D of A. 

In this context, the coproduct comes from the action of a differential operator on a product of two functions. The 
action of di on the product fg is given by the Leibniz rule di(fg) = (dif)g + f(dig). For a product ot two partial 
derivatives we have 

didjifg) = (Oid j f)g + f(a i d j g) + (d i f)(d J g) + (a j f){d i g). (17) 

More generally, for any differential operator P G A, we can write P(fg) as a sum of terms that are the product of 
a differential operator acting on / and a differential operator acting on g. We write this using Sweedler's notation 
P{fg) — S(^ > (i)/)(-P(2)5)- F° r example, if P = di we have a sum of two terms, in the first term P (l) = di and 
P (2) = 1 (with the convention that, for any function /, If = f) and in the second term P (1) = 1 and P (2j = di. The 
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idea of the coproduct is now to remove the reference to the functions / and g and to keep only the sum of terms 
with P (1) on the left and P (2) on the right. This is done formally by defining the coproduct A from A to A <8> A as 
AP = P(i) ® P(2)- From the known properties of the action of a differential operator on a product of two functions 
we deduce the following properties of the coproduct: Al = 1 <g> 1, Adi = di ® 1 + 1 ® di and the recursive relation 

A(PP') = £(PP')ci, ® (PP%> =EE i V> i w® P e> i M- 

From the last rule we obtain A(didj) = {didj) ® 1 + 1 ® (didj) + dj <8> + 3j (g) 9j, and we recover equation (fT7|) . The 
main property of the coproduct is its coassociativity, which means that 

^A(P (1) )®P (2) = ^TP (1) ®A(P (2) ) = ^P (1) ®P (2) ®P (3) . 

For example, if P = 9j, 

p w ® p (2) ® p (3) = ^( A9l )® 1 + ( A1 )®' 91 

= d i <g>l®l + l<g)d i <»l + l<8>l®d i . 

With this definition we can obtain the action of P on a product of three functions as P(fgh) = 5Z(P(i)/)<8>(P(2)flO® 
(/'.. /')• 

After this introduction, we can now define the algebra of functional derivations. The main changes are that the 
partial derivatives are replaced by functional derivatives with respect to external sources, and the fact that the 
anticommutativity of external sources generates signs in the formulas. 



B. The algebra structure of V 

The symbol 8 is used to denote the functional derivative with respect to the external souces rj(x) or fj(x). More 
precisely, since the external sources are two-dimensional vectors, 8 stands for the functional derivative with respect 
to rj s (x) or fj s (x), where s — 1 or s — 2. Products of symbols stands for repeated derivations. For instance, if 
dx = S/Srj 1 (x), 82 = S/dfj2(y) and 83 = 8/Srj2(x), then 

A3 

818283 



Srji (x)5ri 2 {y)5ri2{x) 

The functional derivatives anticommute, thus 88' = —8'8 for any functional derivatives 8 and 8'. Therefore, for any 
functional derivative 8, we have 88 = 0. 

A basis of the vector space D of functional derivatives with respect to external sources is given by the products of 
derivations 8\ . . . 8 n for all n > 1 and the unit 1. Here, the unit is not the constant function 1, it is a symbol that 
satisfies 18 = 81 = 8 for any functional derivative 8. Thus, for instance, 

41 + 2^ + i S - 

Srji(x) 6 Sr] 2 (y)Sr]2(x) 

is an element of T>. 

In T>, the terms of the form d\ . . . d n generate a subspace of T> denoted by T> n (for n > 0). The elements T>o have 
the form Al, where A is a complex number. If D 6 T> belongs to T> n for some n, we say that D is homogeneous and 
its degree, written deg(Z?), is n. For instance deg(l) = 0, deg(<9) = 1, deg(<9<9') = 2. The vector space V becomes an 
algebra if we define the product of two elements of T> to be the composition of derivations. For instance, the product 
of 8\ and 82 is 8182- This product is anticommutative. It can be checked that T> is an associative algebra with unit 
1. Moreover, deg(DD') = dcg(D) + dcg(D') for any homogeneous elements D and D' of T>. From the degree deg(Z?) 
of a homogeneous element D we can define its parity \D\ by \D\ = if deg(D) is even and \D\ = 1 if deg(Z?) is odd. 
If |-D|=0 (resp. \D\ = 1) we say that D is even (resp. odd). 

Now we prove a useful property of the product in T>: if D and D' are elements with a specific parity \D\ and \D'\, 
then 

£>£>' = {-l) mD '\D'D. (18) 

An important consequence of this is the fact that an even element of T> commutes with all elements of T>. To prove 
equation (|18fl . we first show it for homogeneous elements. We start with D = 8 and D' = 8[ . . . d' n , then 88[ . . . d' n = 
(—l) n 8[ . . . 8' n 8 because 8 must jump n times over a 8' . Now, if D = d\ . . . d m) 8 m jumps over D' , giving (—1)™, then 
8 m ~i jumps over D' giving another (-1)", and so on until di and we obtain DD' = (-~l) mn D'D = (_i)deg(£>)de g (£>'). 
Equation l|18|) is recovered because (— l) de s( z> ) de s( £) ) — {—\}\d\\d | jf £) an( j jy are n0 ^ homogeneous but have a 
definite parity, they can be written as sums of homogeneous elements, and the result follows by linearity. 
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C. The coalgebra structure of T> 

We introduce now the coproduct A of T>. In concrete terms, the coproduct of an element D of D is the sum of the 
ways to split D into the product of two elements of V. Formally, the coproduct is defined as a map from V to V <8 V, 
where ® stands for the tensor product. We recall the main property of the tensor product: for any D, D', E,E' G T> 
and A, A', p, fx' e C, 

(\D + X'D') ® (fiE + n'E') = XfiD ®E + Xfi'D <8 E' 

+X'iiD' ®E + X'u'D' (8 E'. 

The coproduct of the elements of smallest degrees is given by 

Al = (19) 
Ad = d®l + l®d. (20) 

To define the coproduct of elements of higher degree, we need a notation for the coproduct. Following Sweedler, we 
write AD — ^2D W ® -D (2) . For instance, if D = d, the sum has two terms. The first term is D {1) = d, D {2) = 1 the 
second term is D (l) = 1, D {1) = d. The coproduct can now be defined recursively by 

A(DD') = ^(-^'"^"^'(^(x)^))®^)^))- (21) 

As an exercise, we calculate A(dd'), so that D = d and D' = & , Equation QQJ gives us Ad = 9(81 + 1(89 and 
Ad' = & (8 1 + 1 <8 3'. The first term of A(99') is obtained from formula with = 9, D (2) = 1, D' (1) = d' and 
Z?' 2) = 1. The degrees are \D (2) \ = 0, \D' (1) \ = 1 and their product is |£> (2 )||D' 1) | = so we obtain the term dd 1 <& 1. 
The other terms are calculated analogously and the result is 

A(dd') = dd' ® 1 + 1® dd' + d® d' - d' ® d. 

The minus sign is due to the fact that the corresponding term comes from D (1) = 1, D (2) = d, D' (1) = d' and D' (2) = 1, 
so that |D (a) 11.0^,1 = 1. 

It can be checked^ that the coproduct of a basis element D = d\ . . . d n of V is 

AD = D® l+l® D 

n-l 

+ ^2 X/( _1 -' Cr ^ <T ( 1 ) • ' ' ^ ct (p) ® ^(P+1) ■ ■ ■ ^(n)j 

p— 1 cr 

where a runs over the (p, n— p)-shuffies and (— 1) CT is the signature of the permutation a. Recall that a (p, n— p)-shuffle 
is a permutation cr of {1, . . . , n] such that <r(l) < cr(2) < • • • < cr(p) and cr(p+ 1) < • • • < cr(n). Notice that we always 
have D = D {1) D {2) . 

With this definition, we know the coproduct for a basis of T>, the coproduct of a general term of V is obtained by 
linearity: A(AL> + X'D') = X(AD) + A'(AD')- 

The most important property of the coproduct is its coassociativity. We saw that the coproduct of an element D 
gives the ways to split D into two elements D {1} and D m - Now assume that we want to split D into three elements. 
We can achieve this either by splitting D m or by splitting D {2) . Coassociativity means that the result does not depend 
on this choice. This is expressed more formally by (Id <S> A) A = (A (8 Id)A. For example the reader can check that 

(Id<8A)Al = 1 (8 1 <8 1 = (A (8 Id)Al, 

(Id (8 A)A<9 = 9(81(81 + 1(85(81 + 1(81(85 = (A <8 Id)Ad. 

The coproduct A is coassociative for all elements of T> 57 . It can also be shown that the coproduct satisfies AD = 
X)-O(i) <8 D {2) — ^2(— 1)' D < 1 )"- D < 2 )'l? (2) (8 D m (this property is called graded cocommutativity). 

We can define recursively the splitting of D into n parts by A (0) D = 1, A^D = D, A^D = AD and A^D = 
(A (8 ld n ~ 2 )A^ n ~ r W for n > 2. The result of the action of A^ 1 ) on D is denoted by 

AWfl = D m (8 •••®-D (n) . (22) 

To make a Hopf algebra, we need also a counit and an antipode, but we shall not use these concepts in the present 
paper. 
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D. The derivative of a product 

To show immediately the power of the Hopf algebraic concepts, we prove the following formula for the derivative of 
a product of two functions. If D G D is a product of functional derivatives and u and v are functions of Dirac fields 
and sources we have 

D{uv) = J2(-^ Di2)M (D w u)(D l2) v). (23) 

In this equation, \u\ is the parity of the function u. The parity of a function is defined as follows. We first define the 
degree of a function: for a Dirac field or a fermion source we have deg(V>) = deg(V0 = deg(?7) = deg(?7) = 1. The 
degree of a product of fields and sources is the sum of the degrees of the fields and sources: deg(ui>) = deg(u) +deg(w), 
and the parity of a function of fields and sources is equal to the or 1 when its degree is even or odd. Notice that, 
if deg(D) < deg(u) we have \Du\ = \u\ + \D\ modulo 2 because deg(Du) = deg(u) — deg(D). The proof of ffify is 
recursive. Equation H23|) is true for D = 1 because l(uv) = uv and for D = d because of Leibniz' rule JHJl. If this is 
true for all elements of degree up to n, take D an element of degree n and define D' — dD. On the one hand 

D'{uv) = d{D{uv))=Y J {-^ D ^ M d({D {1) u){D {2) v)) 

= J2(-l) lD ^ lul (dD (1)U )(D, 2) v) 

+(-l)l D wlW+l^wl+l«l(D (I) «)(flD (a) «)). (24) 

To obtain the last line, we used Leibniz' rule and the fact that |-D(i)it| = |u| + \D^\ modulo 2. On the other hand, 
by equation (|2l"1) 

A(clD) = ^(3D (1) )®D (J) + (-l)l D (»l]) (1) ®(9]) (2) ), 

So that, if equation b'2'Mi is true, 

D'(uv) = J2(- l ) lD(2>lH (9D {1) u)(D (2) v) 

+ (-l)\ D w\+( 1 +\ D (2i\)H( D(1)U )(dD (2) v). 

But this is indeed equal to (J2U, so equation (|2*3^l is satisfied for D'. Since the elements dD generate 2?„+i, equation 
is true for V. 
More generally 

D( Ul ...u n ) = ^(-l^^Ef^l^wlKI^^)...^^). (25) 
The recursive proof is left to the reader. 

E. Elimination of closed shells 

As a second application, we calculate tr(pN ) when the system is composed of closed shells and open shells. A 
closed shell is an electron state ik which is occupied in all states \K). The open shells are the electron states which 
are present in some but not all states \K). Thus, the closed and open shells have no electron state in common. We 
rewrite equation l|16|) as ti(pN°) = e d {uv) where d = J2 n d 2 /da n da n , u = a mi a mi . . . a mc ,a mc , describes the closed 
shells containing C electrons and v — J2kl PLK&jt&i! ■ ■ ■ &j N &iN describes the open shells. Notice that in u the index 
of each a is the same as the index of the following a. This is because the electron states are ordered so that the closed 
shell have an index smaller than the open shells, and the closed shells are occupied in all \K) and \L). To calculate 
tr(/5A°) we first compute d(uv). According to equation llL'Mli 

d(uv) = ^(-l)l^)H"l(d (1)M )(d (2) «)=^(d (1)U )(d (2) «), (26) 

because |u| = 2C = modulo 2. Now 
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The terms du/ da n dv / da n and du/da n dv/da n in equation l|26[l are zero because the closed and open shells have no 
state in common. Therefore d{uv) — (du)v + u(dv). Moreover, 

c 

du = > 



fc=l am k a "i k 



is a sum of closed shells, so we can apply the same argument again to show that 

1=0 ^ ' 



Therefore 



oo oo k 1 



kl y ' ^^l\(k-l)\ 

k=0 k=Q l=0 K ' 



11 

T\ dlu E — dmy = (^X^) = U( l + a mi a mi )(e d v). 



1=0 m=0 i=l 

In other words, the closed shell factorize in ti(pN°). This result will be important to restrict the size of the problem. 

Notice that, in the proof, we used only the fact that the closed and open shells have no electron state in common. So 
the same reasoning shows that, if the system is composed of two independent subsystems, then N^ L is the product of 
the N^ L of both systems. More precisely, if all states can be written as \K) — \K\) A \K%), where A antisymmetrizes 
the electron states of \K\) and I.K2), where \K{) has the same number of electron states for all \K) and where no \K\) 
and |if 2 ) have any electron state in common for any \K) and \K'), then Nf CL — N^ Ll N^ L2 , 

VI. CALCULATION OF W° 



It will be very useful to define W° = log(Z°). If the system has N + C electrons with C electrons in closed shells, 
Zp can be written 

c 

Zp = exp[-i / fj(x)G^(x,y)r](y)'^dy]\J_(l + a mi a mi )p(a,a), (27) 

with 

JV 

p(a,a) = ^2p k (a,a), (28) 

k=0 

where pk contains products of k a and k a. More explicitly 

p N (a,a) = '^2pj N ...j 1 ,i N ...i 1 a jl a il . . .a ]N a lN , (29) 

(N — k)l V ^ da n da n . 



1 (sr^ a \ N ~ k 



In particular, 



po(a,a) = tr(p), 



px{a,a) = ^Pj N ...j 1 ,iN--A 1 yJ2ot jl a il Y[Sj p ,i p 

1=1 pjti 

JV-1 N 

+e e (-i) i+m ^«^n^^ n ik 

1 = 1 m=l+l p<l Kp<m p>m 

N l-l 

+EE(- 1 ) !+ % q -I1u n sj^Us^), (si) 

1=2 m=l p<m m<p<l p>l 
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will be useful. It is important to isolate pi(ot, a), which depends linearly on a and a, because it will become a part 
of the free propagator. 

The closed shells are dealt with easily: 



log ( JJ(1 + )) = ^log(l + a mi a mi ) 



1=1 

C oo 



t ^ ( — 1)™ +1 

2=1 n=l 

However, a mi and a mi are fermionic variables, thus (a mi ot mi ) 2 = a mi a mi a mi a m; = — a mi a mi a mi a mi — because, 
as fermionic variables, <5^. = a^. = 0. Consequently, only the term n = 1 remains in the sum and 

c c 

log (n^ 1 + 

i=l i=l 

This result is important because it justifies the fact that the propagator of the Green function in many-body theory 
is obtained by summing the contribution of all occupied shells. We see now that this procedure is justified when the 
vacuum |$o) can be written as a full shell. In all other cases, this procedure must be modified. The modification 
comes from the term p(a, a) that we write 

N N n. 

k=l k=l T \P' 

The usual convention is to impose tr(p) = 1, but we want to relax this constraint for later convenience. Thus 

N (- \ 

log(p(a,a)) = log(tr(p))+log(l + ^M 7 ^) ) 

= log(tr(p)) + ^i^ + p c (a,a), 
tr(p) 

where p c (a, a) is defined by the last equation. We can write p°(a, a) as 

oo 

p c (a,a) = ^2p c n (a,a), (32) 

n=2 

where p c n is the sum of the terms of p c which have degree n in a and degree n in a. Notice that the sum over n is finite. 
For instance, if the states \K) are built by choosing TV electron orbitals among M (for instance, for Cr 3+ , we have 
three d electrons so that N = 3 and M = 10). Therefore, p n (a, a) 1+M ^ n = because in each term of p n (a, a) 1+AI ^ n 
at least one a Ip is found twice and af p = 0. Therefore, /0^(<5, a) — for n > M. 
If we gather all these results we obtain that 

c 

W° = \og{Z° p ) = -i / fj(x)G%(x,y)r)(y)dxdy + ^a ro< a TOi 
J i=i 

+ log(tr( / 5)) + ^i^+p c (a,a), (33) 
tr(p) 

where we recall that a n = J(fj + (x) — fj^(x))u n (x)dx and a n — J u n (x){rj + (x) — r]-(x))dx. The term containing 
Gq(x, y) is linear in fj and 77. Thus, we shall include the other linear terms by defining 

C , s , 

G° p (x,y) = G*(x,y)+i{Y,u mi {x)u mi (y) ' ^"""'^ 1 ' _l 



i=i 



tr(p) ; V - 1 1 



with pi{x,y) defined so that pi(a,a) = J Tid{x)pi(x,y)rjd{y)dxdy, in other words, p\{x,y) is obtained by replacing 
all a.i k by Ui k {x) and all ctj k by Uj k (y) in equation (|31|) . It is at this stage that, when the system has only closed 
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shells, the effect of the closed shells is entirely taken into account by adding the occupied orbitals to the free Green 
function. This procedure, which is universally used in the quantum many-body approach, is usually deduced fromt 
the particle-hole transformation. This transformation is itself justified by showing that the Hamiltonian without 
interaction Ho is left invariant (up to a pure number) 32 . However, this justification falls short of being a proof that 
this procedure is valid at all orders of the interacting theory. From the previous discussion, we see that the procedure 
is correct at all orders when the noninteracting system can be described by a single Slater determinant (i.e. a closed 
shell). However, the most interesting phenomenon occurs when open shells are present. We rewrite 

W° p = -if rj(x)G° p (x,y)ri(y)dxdy + log(tr(p)) + p c (a,a). (34) 

This is the final result of the section. 



VII. THE GREEN FUNCTION HIERARCHY 



In this section, the Green function hierarchy is established in the presence of open shells. 



A. Definition of Green functions 

According to the discussion of III Al the expectation value of the Heisenberg field iPh(x) is given by 

{^h{x)) p = _i s ^L\ \v=v=o- ( 35 ) 
The density matrix is normalised by tr(p) = 1, so that Z p \^ =ri= § = 1. Therefore, we can also define 

= (iiSy)'-'<- (36) 

Although these definitions are equivalent, equation 13(j|) has some advantages over equation (I35f) : (i) If we multiply 
p by A, equation (|3^|) is not changed because the factor A is cancelled between the numerator and the denominator. 
Thus, it is possible to relax the constraint tr(p) = 1 and we are enabled to consider unconstrained density matrix. 
In particular, we can use p = exp[— f3H] for equilibrium quantum field theory, (ii) If equations (|35|l and l|3tj[) are 
written as a sum of Feynman diagrams, equation (|35|l has vacuum diagrams which are cancelled by the denominator 
of equation i|36|) . in other words, only equation i|36|) is a sum of connected diagrams, (iii) When the density matrix 
p is that of the vacuum (i.e. p = |0)(0|), equation (|36[) is the Gell-Mann and Low equation 3 ?- which is known to be 
correct, (iv) Equation l|36ll has been used successfully since the early days of nonequilibrium quantum field theory^. 

It turns out that a complete set of equations cannot be obtained by functional derivatives with respect to rj + and 
fj+ alone. So we define the following expectations values: 



GM*)>, = (feW>P = (^)l«=o> 



B. Hierarchy of disconnected Green functions 



We rewrite equationlO as Z p = c lD Z° 



p 
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The operator D contains products of 2 or 4 functional derivatives, thus D is even and D commutes with the elements 
of V. Thus, if /3 = r]±(x) or (3 = fj±(x) 

5(3 5(3 ' 



We use the fact that Z° = e w ° with \W°\ = to get 



5Z P _, n fSW* w o A _ in/ 5W° 



n=0 ' M 

The action of the operator D n is expanded with equation 123|l . using \5Wp/5(3\ = 1: 

5Zp y> H) ra \^f_n]ggolf£)» ^ 
5/3 ^ n! ^ l J 1 (1) (5/3 ^f^' 

n— 

We transform this infinite sum into a finite sum by using reduced coproducts. The reduced coproduct with respect 
to D is denoted by A'D. It is defined as follows, the reduced coproduct with respect to D of D itself is defined by 
A'D = AD -1&D-D&1. The Sweedler notation for it is A'D = J2 D (l/) ® D (2 / y The reduced coproduct of D n 
is defined recursively by 

A'(D n+1 ) = Y,{-^ D ^ DW ^ D ^®^ D v'y (39) 

This is extended to n — by A'(D°) = 1 ® 1. An equivalent definition is that A'(D n ) is the sum of all terms of 
A(D n ) which do not contain any D. The relation between A(D n ) and A'(D n ) is given by 

n n — k | 

AP") = EE kW(n n J k ^ lV D ^ lDk ® (40) 
This can be shown by a recursive proof. The definition of A'D gives us 

A£> = £)®l + l<g>.D + E i V)® i V)' (41) 

so equation (|40|l is true for n = 1. Assume that it is true for all D k for all k up to n. From equations (J2TJ, l|4U|) and 
(gTJ we obtain (using |£>| =0), 

A(D^D) = E 

A.'+/+m— n 

+D™, } D k ® £>^)£) z+1 

+ (_l)l D d')ll^')l£)^ )J D (1 , )J D fe (g, J D^ ) £» (2 , )J D z ). 

Using the recursive definition (|39() we get 



+D™, } D k ® D™ 2l) D l+1 + D^ l + 1 D k <g> L*™,} 1 ^ 



This can be rewritten 



v ' ^ \(k-l)\l\m\ kl(l-l)lml k\ll(m-l)\J 

k+l+m=n+l y ' y ' y ' 

D™ l} D k ®D™ t) D l . 
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The first three integers can be summed to 

A(D«D)= £ ^^^^^ ! , 

fc+Z+m— n+1 

and equation (|40|) is proved for _D n+1 . 

By summing equation (|4U|1 over n we obtain the important identity 

oo 1 

Ae D = ^- J D^ ) e D ® J D^ ) e D = (A'e D )(e D ®e 15 ). (42) 

n=0 

Note that this identity is true for any graded commutative Hopf algebra and any D of degree > 0. 
Using identity l|42|). the equation 138fl for 8Z p /8(3 becomes 

w = E^r : E(- 1 ) |Dr2 ' )l (^,^^)W)^^ P ), 

H 71=0 ' ' 

where = e~ lD W p adds the electron-electron interactions to the cumulant W° of the moment generating function 
Zp. Since the cumulant W p is a finite polynomial in a and a, the interacting cumulant W p is also a finite polynomial 
in a and a. Now each D™, (for contains at least to functional derivatives with respect to fj± or rj± (this is 

why the reduced coproduct was defined), thus DVf,. is zero for to large enough. In fact, to = 2M — 1 is a possible 
bound and we obtain our final formula, isolating the contribution of n = 

(43) 

We have transformed the infinite sum l|38|) into the finite sum (|43|l . To be complete, we still have to replace the 
disconnected Green functions defined by functional derivatives with respect to Z p by connected Green functions 
defined by functional derivatives with respect to W p = \ogZ p . 



C. Calculation of W p 

Apparently, W p = e~ lD W p includes some interaction in W p , but in the interaction Hamiltonian H lnt that we 
consider, we have W p = W p . Indeed, these contain integrals over d = S 2 / 5rj±(x)Sfj±(x) . The action of d on the term 
containing the Green function G p (x,y) is irrelevant because it gives a term independent of a and a. For the action 
on p c (a,a) we have 

Sp c _ dp c 5a n _ ^ dp c ^ 
8fj±(x) t-f da n 5fj±(x) ^ da n 



A 2 n c fl 2 n c 

P _ ST^ ° P - / \ u / \ 

5r]±(x)6fj±(x) da m da„ 



Remark that the right hand side of the last equation does not depend on the sign ± of the source. The diffential 
operator D can be written as D = D + — _D_, where D+ and D_ are the same operators, but the first one involves 
derivatives with respect to the + sources and the second one with respect to the — sources. According to our remark, 
D+p c = D-p c . Thus, Dp c = and W p = W°. 
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D. Hierarchy of connected Green functions 

In formula l|43|) . the differential operator D 1 !,. acts on Z p = e ' = X^^Lo W™/n\. Thus, we must determine the 
action of a differential operator on W™. Notice that \Z p \ = 0, thus \W P \ = 0. 

So we take an even element u (even means that \u\ = 0) and a differential operator d such that deg(d) > and 
we want to calculate du 11 . We shall use now the standard reduced coproduct A defined, for any element d € V by 
Ad = Ad — d<E>l — 1 ® d, and we write Ac? = ^ d a) g> d {1) . This reduced coproduct is coassociative. The basic identity 
that we need is 

d ^ n ) = E (D un ' k E d ^ u ■ ■ ■ d ^ ( 44 ) 



k=l 



where d {1) u = du if k = 1 and A^d = Yl d <X) ® "" ® d (k) is defined recursively from A as in equation l|22|) . For 
example, using equation (|23|) and |u| = 

d(u 2 ) = '^(d {1) u)(d (2) u) = (du)u + u(du) + 22(d(r>u)(d m u), 

— 2u du + y^(d ( i)M)((i ( 2)tt), 

and equation (|44|) is valid for n — 2. The general case is proved recursively. Assume that it is true up to n, then 
d{u n+1 ) = ^2(d a) u n )(d i2) u) = u n du + d{u n )u + ^2(d (L) u n ){d {1) u), 

= u n du + E ( T ) un ~ k+1 E ■ • ■ d ^ u 
k=l ^ ' 

+ E ( I E • • • 



fc=l 



i ^ + E (fc) M "~ fc+1 E d <i> u ■ • • <w 
Eli.) u "~ fc E d <i' M • ■ • d (±±i) M 



fc=i 



fc=i 

n+l 



= u"o;m + y \u n k+1 ^2d (L) u...d (k) u 



k=1 
n+l 



£U_i)« B "* fl E d <» u "- d (u u 



= Ert> n+1 "E d a,".-.^,- 



We used the coassociativity of the reduced coproduct. From equation (|44|l we can calculate 

oo oo n j \ 

d ^ = E^") = E^E C j^Ev---^ 



nl A — ' n\ 

n— n—1 k—1 v 7 

= E ^E^E d a> u ---^> M ' 

m=0 ' fc=l 

= eU E w E d <i) M --- d <^ u - ( 45 ) 

k=l 

The sum over k is not infinite because A^ k 'd = if fc > deg(d) and the sum stops at k = deg(d). More generally, for 
an analytic function f(z), 

d (f( u )) = E fci E d <i) u --- d ^) u ' 
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where f^ k \u) is the fc-th derivative of / at u. The cocommutativity of the coproduct ensures that the factor 
disappears from the expanded formulas. 

If equation 1)45(1 is applied to u = W p , we obtain a relation between unconnected Green functions {\/Z p )dZ p 
and connected Green functions dW p . For instance, if d = dd', then Ad = d <E> d' — d' <S) d and (\/Z p )dZ p = 

dW p + (X/2)(dW p )(d'W p ) - (1/2) (d'Wp)(8W p ). At fj± = r} ± = we obtain dZ p = dW p . Similarly, if d = dd'dd', 

where d and d' are derivative with respect to rj and d and d with respect to fj, we find at fj± = rj± — 0, dZ p = 

dW p - (ddW p )(d'd'Wp) + (dd'Wp)(d'dW p ). 

Equation 145|l is now introduced into 143|) . where we use the fact that = \D n \ — l-D™/,! + l-Dfy-Ji so that 
D n I — I D" I • 

n— 

oo ^ 

xZ "EfciE (^r 2 ',<D^p ■ ■ ■ ) (^, tS ) Wp) ■ 
fc=i 

Using again the definition of W p in terms of Z p , we obtain an equation involving only the connected Green functions: 

dW p 1 dZ p 

~Jf3~ ~ Yplp 

2M-1 , .>„ oo , u\ D 7i')\ AW 1 

n\ fc! V (5/3 ; 

n=0 k=l ' 

(D^Wp)...(Dl, m W p ). (46) 
This sum is finite because, for each n, the sum over fc stops at fc = deg(D™, ) ). 



VIII. CONCLUSION 



This paper had two purposes: (i) to determine the hierarchy of Green functions for degenerate systems, and more 
generally for systems whose initial state cannot be written as a Slater determinant; (ii) to show the power of quantum 
groups and Hopf algebras to solve problems of quantum field theory. A detailed application of the formulas obtained 
in this paper can be found in reference^. 

In this paper we dealt with a nonrelativistic electronic system with Coulomb interaction. A generalization to 
quantum electrodyamics is possible, which would provide an alternative to the new methods recently developed to 
carry out quantum electrodynamical calculations of many-electron systems£2i§2i£ii§2iS. Again, the present method 
has the advantage of being self-consistent and of preserving the symmetry of the system. 

Moreover, a functional derivation of the energy with respect to the density matrix provides equations that enable us 
to unify the Green-function formalism and the diagonalization method of many-body theory. This will be presented 
in a forthcoming publication. 
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X. APPENDIX : CALCULATION OF THE TRACE 



The calculation of tr[pN Q (rjd, rid)] is an essential ingredient of this work. We rewrite the density matrix as p = 
^2kl Plk\L){K\, where \K) and \L) are Slater determinants defined by \K) — b\ N . . -frjjO) and \L) — bj N . . .frjjO). 
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Here b\ k and b'^ are creation operators of the one-electron orbitals indexed by ik and j\. The indices are ordered 
(«!<•••< i jy ) ii < ■ ■ • < JaO- The total number of electrons in the system is N. Moreover, |0) is the true vacuum 
of the system (i.e. containing no electron). We must calculate trfp./V ^, rjd)] — J2kl PlkN^ l with 



N° KL = (Jf|:exp(i / fj d (x)1>{x)+${x)ri d {x)dx):\L). (47) 
The fields are expanded over time-dependent eigenstates of the one-body Hamiltonian 

n n 

where u n (x) are the time-dependent solutions defined in section TlV Bl and n is the index of the electron orbital, 6 n ,^n 
arc the annihilation and creation operators of an electron in orbital r>S^. 



We can rewrite N^ L as 



00 'l r r 1 

V ' : ' = EtT^I : (E / nd(x)u n {x)&xb n + bi / u(x)r] d (x)dx) :\L) 



"KL 

/-II 



i=0 



where a n = J rj d (x)u n (x)dx and a n = / u n (x)r)d(x)dx are anticommuting variables. To calculate N^ Ll we first notice 
that the anticommutativity of b n , W n , a n and a n for the normal product gives us the commutation rules :ctibi(Xjbf = 
■.ctjbjCiibi:, ifijbjfctay: = MajCiibi: and :b\a.ib^otj: = :b^otjb\ai:. Thus, we can expand the power with the binomial 



formula 

00 1 1 



1=0 k=0 x ' 

E ( K \ b li a m ■ ■ ■ b l k a n k a nk+1 b nk+1 . . . a ni b ni \L) 



m— ri| 

00 ./ / 

— i l 



*-n k <~*-n k + 1 



= E^EU (-D fe+i( '- 1)/2 E 

(K\bl...bl k b nk+1 ...b ni \L). 

The transition between \K) and |L) is zero if I ^ 2k or if I > 2N because \K) and \L) contain N electrons. Thus we 
obtain the finite sum 

" (-1)* 



n kl - E-pp- E 



fc=0 



(if|6t 1 ...6j lk 6 mi ...6 ro jL). (48) 
A. Hopf calculation 

Hopf algebraic techniques will be used to obtain an explicit expression for N^ L . We first denote 

Arl = {K\tf ni . . . b\ k b mi . . . b mk \L), 

and we write u = bi x . . . bi N , v = bj N • • • &] t , s = b\ ti . . . b\ lk and t = b mi . . . b mk . Thus Akl = (0\u(:st:)v\0) and we use 
the Hopf version of Wick's theorem^! 

(:st:)v = ^(-l)l^)KI^)l + l^)l)(: S t: (1 > (1) ):(: S t: (2) z; (2) ):, 
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Therefore 



In general 



Akl = (u\(:st:)v) , 

(:c 1 ...c m :\:d 1 ...d n :) = 5 m ,„(-l) n(n - 1)/2 det(M), (49) 



where Ci and are creation or annihilation operators and M is the n x n matrix with elements My = (ci\dj)r*. The 
Laplace pairing (ct\dj) is obtained from = <5$ i3 -, (i>i|b 3 ) = 0, (&-|b 3 ) = and = 0. Because of the value of 

(pi\dj), (:ci . . . c n :\:d\ . . . d n :) is zero if any c; is a creation operator or any dj an annihilation operator (because one 
row or one column of M is zero). Therefore, we need s (1) = 1 and i (2) = 1, so that s (2) = s and t (1) = t: 

Akl = ^(-l^^'^'l^h^H:^:), (50) 
= ^(-I)l-< 1 )IN+I t ll s l+I"( 2 )IN(^ (1) )( U(1) | S )( U(2) | W(2) ). (51) 

We rewrite v = (-l)-^" 1 )/ 2 ^ ...b] so that 



31 3n 
N 

Aw = ^^(-ir6 Ml) ...6 Mp) ®6 Mp+1) ...6 MA . 



P=0 CT 

N 



Av = (-i) Ar ( jv - i )/ 2 yy(-i)^ ... 6 t 86 t ... 6 t 

q=0 r 

where er runs over the (p, AT — p)-shufnes and r over the (g, A~ — q)-shuffles. A (p, N — p)-shufflc is a permutation a 
of (1, ... , AT) such that <r(l) < a (2) < ■■■ < a(p) and er(p + 1) < • ■ • < a(N). If p = or p = N, a is the identity 
permutation. Equation l|49|l applied to (|51|l . gives us p = k and = fe so that = |s| = \t\ — k, \u (2) \ — N — k and 

A KL — 'y^^_^N(N-l)/2+(N-k)k+k(k~l) + (N-k){N-k-l)/2 

^(-l)^det(<5 mpjT(g) )det(^ (p) ,„Jdet(^ (p)jT(9) ), (52) 

(TT 

where p and q run from 1 to A; in the first two matrices and from k + 1 to A" in the last one. The determinant of 
a n x n matrix ay is det(a) = X)a(~1) A T[i=i a iX(i) — X)a(~1) A Yl7=i a Mi )j> wnere A runs over the permutations 
of n elements. Therefore, to calculate the last determinant in equation (|52|l . we must sum over all permutations of 
r(k + 1), . . . , t(N), but the indices satisfy i\ < ■ ■ • < ijv and ji < • • • < jw- By definition of the (fc, A" — /c)-shuffle, 
we have V(k+i) < • • • < 2 CT (Ar) and j T (fe+i) < • • • < i T (jv) so any permutation of r(fe + 1),. . . ,t(N) would break this 
ordering (for example, $j u i a $j 3 ,ii — because j\ < j% and i\ < Thus the only nonzero term of det(#j , p 5 , 3 , , ) is 
^ IH1) ,jV (H1) ■ ■ -<> MN )j T (iv r This gives us the following expression for A KL : 

Akl = ^(-l) fc(fc - 1)/2 ^(-l) CT+r det( ( 5 mpjT(9) )det(^ (p) ,„J 

err 

n w<„. 

where a and t run over the (k, N — k) shuffles. 
To calculate det(S mpt j T , ,) we write 

A 

where A runs over the permutations of {1, . . . , k} and we obtain 

a TOl . . . a mk det(<5 mp Jt(?) ) = fc!<x, T(1) . . . <x, T(fc) . 

m 1 ,...,m k 
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Hence 

JV 

N° KL = E(- 1 ) fe(fe+1)/2 E(- 1 ) ff+T «M 1 ,---«v (fc) 

k=0 err 

«Jt W • • • »3t(«0 det ( 5 V(p)Jr W )- 

Therefore, our final result is 

iV fe fe JV 

^ = E(- i ) fc( ^ i)/2 E(- i ) CT+T n^ (P )ii«M P , n 

JV fe JV 

= EEMr^n^,^,) II ^ w jr W . (53) 

k=Q or p=l p=k+l 

where we recall that a and r run over the (fc, N — k) shuffles. 

We calculated N^ L for a system where all the states have the same number of electrons, but the same methods 
can be used when \K) and \L) have a different number of electrons. 



B. Alternative formula for N% L 

Now we are going to derive the alternative formula for N^ L . This result can be obtained directly from 
eauation (|53|l but we shall provide an independent proof. 
We first rewrite the expression (|48|) for N^- L as 

N (-l) k 

N °kl = ^ { j^{0\b il ...b iN B k + BHl...b] i \0), 

where B — Y^n^nbn an d B + = ^2 n b\fit n . It is easy to prove recursively that \b iN ,B+] = ka iN B k ^ 1 and [B k bj N ] = 
kB k ~ 1 ctj N . If we write \K~) = b\ . . .6^, so that (K\ = (K~\bi N , we obtain the recursion 

{K\B%B k \L) = (K-\B k ;b lN B k \L) + k(K-\a iN B k ;- 1 B k \L). 

If we use now \L~) = b^ Ni . . . bt, we obtain the following recursive equation between the matrix elements of B+B k 
for TV particles and N — 1 particles: 

(K\B k + B k \L) = (K-\B k B k \L-)S jNtiN 

-h?(K-\BX- 1 B k - 1 \L-)a js a iN 

-k{-l) N -HK-\Bl- l b] N B k \L-)a lN 

+k(-l) N - 1 (K~\B k + b tN B k - 1 \L-)a jN 

-{K-\B%b\ N b iN B k - x \L-). (54) 

Now we are going to show that the expression (jlHfl satisfies the same recursive equation. We write d = ^2 n d 2 /da n da n , 
u = (cti 1 aj 1 . . ■ ai N _ 1 aj N _ 1 ) and v — cti N aj N , so that N^ L — e d (uv). Equations l|23|) and (|42ll yield 



oc 



v 

p=0 y 



The sum is not infinite because e d v = Sj N! i N + ctj N a.i N so the sum stops at p = 2. Using 

A'd — ® ® ® ® ® ® 

' da n da n da n da r , 



2 Q2 q2 



A'd 2 = V 9 & <>2 ° 



da m da n da m da n da m da n da m da n 

Q2 Q2 Q2 Q2 

+ ^ 7^r- ® — ^ b 



da m da n da m da n da m da n damdan' 
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we obtain the recursion 



dr \ / rf \ r t d \ — d(e d u) d(e d u) 
e (uv) = (e u)5 jN ,i N + (e u)a ijv a iiv + — a iN + — — 

d 2 e d u 



da JN N da iN 



a 3N 



We use the derivatives 



to obtain 



da lN da jN 



9a JN JJV da, 

d{K-\B\B k \L-) 



da jN 



W-ikiK-lBl-H^&lL 



Similarly 



d{K-\B%B k \L-) 
da iN 

d 2 {K-\B% +1 B k+1 \L-) 



(55) 



9 * iN 8a jN (k + lf { K-\Bl b l biN B«\L-), 

With these identities, it is easy to show that N < L_ L _ and N^ L satisfy the same recursion as e d (u) and e d (uv). When 
there is only one electron (N = 1) it is easy to show that Nj^ L = Sj ly i 1 + aj 1 ai 1 = e (a^a^). Thus we have 
n kl = e d (a 3l a H . . . a 3N a lN ) for all N. 

Notice that equation l|55|l enables us to calculate N^ L for a closed shell, i.e. when N = M. In that case there is only 
one state \K), where all orbitals are filled, and we have if. = jk for all k = 1, . . . , N. Since all the orbitals are different, 
e d u does not contain im = Jn- Thus the partial derivatives are zero and we obtain e d (uv) = (e d u)(l + cti N ai N ). For 
JV = lwe have N^ L = 1 + cti 1 a il , thus 

N 

Kl = a lk a lk ). (56) 

k=l 
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